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Abstract 

Many new physics models predict sizable modifications to the SM Zbb couplings, 
while the corresponding measurements at LEP and SLC exhibit some discrepancy 
with the SM predictions. After updating the current results on the Zbh coupling 
constraints from global fits, we list the observables that are most important for 
improving the Zbb coupling constraints and estimate the expected precision reach 
of three proposed future e"''e“ colliders, CEPC, ILC and FCC-ee. We consider both 
the case that the results are SM-like and the one that the Zbb couplings deviate 
significantly from the SM predictions. We show that, if we assume the value of 
the Zbb couplings to be within 68% CL of the current measurements, any one of 
the three colliders will be able to rule out the SM with more than 99.9999% CL 
(5cr). We study the implications of the improved Zbb coupling constraints on new 
physics models, and point out their complementarity with the constraints from the 
direct search of new physics particles at the LHC, as well as with Higgs precision 
measurements. Our results provide a further motivation for the construction of 
future e"''e“ colliders. 
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1 Introduction 


The LHC has just started to run at an unprecedented center-of-mass energy, 13 TeV, and 
will be able to probe new physics (NP) at higher energies. At the same time, precision 
measurements of electroweak physics at future e’''e“ machines will also offer powerful 
probes of Beyond the Standard Model (BSM) physics. 


The next lepton collider may not be too distant in the future from us. Several com¬ 
pelling plans exist, including the International Linear Collider (ILC) [^, FCC-ee, formerly 
known as TLEP [^, and the Circular Electron-Position Collider (CEPC) [^. With the 
discovery of the Higgs boson, the primary goal of such future e’''e“ colliders will be to 
produce a large sample of Higgs boson events at around ~ 250 GeV to precisely measure 
the Higgs boson’s properties, acting as a “Higgs factory”. On the other hand, such e+e“ 
colliders could also collect large amount of data around the Z-pole, producing several or¬ 
ders of magnitude more Z-bosons than what was produced at LEP. The large amount of 
Z-pole data would greatly improve the measurement of the several electroweak precision 
observables (EWPOs), which could provide strong constraints on NP. 


Several studies on the measurement of EWPOs at future e'''e colliders have been 
performed @ 0 . mainly focusing on the oblique corrections parameterized by the Peskin- 
Takeuchi parameters S and T [^. NP could also have sizable non-universal corrections. 
The corrections to the Zbb vertex is particularly interesting, and quite generic in NP 
models. For example, being the left-handed top and bottom quarks in the same electro¬ 
weak (EW) doublet, new physics that couples to the top quark usually also affects the 
Zbb couplings [^. Composite Higgs models with light top partners usually predict a large 
correction to the Zb]J)L coupling, unless it is protected by some symmetry analogous to 
the custodial symmetry that protects the weak isospin [^. Additionally, heavy Higgs 
bosons typically couple more strongly to the third generation quarks and modify the Zbb 


coupling through loops 11 


The story is even more interesting on the experimental side. At LEP, the left and 
right handed Zbb couplings are mainly determined by two measurements at the Z-pole: 
R^, the ratio of the Z ^ bb partial width to the inclusive hadronic width, and Apg, the 
forward-backward asymmetry of the bottom quark. The measured value of agrees 
with the most recent two-loop calculation of its SM prediction within 1 a 
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instead, exhibits a long-standing discrepancy with the SM prediction with a significance at 

In addition, SLD directly measured the bottom quark asymmetry 


around 2.5(7 12,14 
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measured value 

SM prediction 


0.21629 ± 0.00066 

0.21578 ±0.00011 

.0,6 

^FB 

0.0992 ±0.0016 

0.1032 ±0.0004 

Ab 

0.923 ± 0.020 

0.93463 ± 0.00004 


Table 1: The measured values and SM predictions of R^, Apg and At according to the 
most recent result from the Ghtter group 14 


with longitudinal beam polarizations, At, which is consistent with the SM prediction 
within 1 a but slightly prefers a shift in the same direction as the LEP Apg measurement 
does. The measured values and SM predictions of Apg and Ab are summarized in 
Table 1 b To obtain the desired modihcation for Apg without violating the experimental 
constraint on and At, a simultaneous modihcation of both the Zbibi and Zb^bji 
couplings is required. To obtain the best estimation for the preferred values of the Zb]J)L 


and Zbjibii couplings, a global ht to all precision data has to be performed (see e.g. 16 17 
for earlier studies). 


A future collider offers great opportunities for further studies on the Zbh cou¬ 
plings. With the huge improvement on statistics at Z-pole, it will surely have the potential 
to resolve the Apg discrepancy at LEP. If the result agrees with the SM predictions, a 
e’''e“ collider can provide very strong constraint on NP models; if the LEP Apg discrep¬ 
ancy does come from NP, a collider will have the potential to rule out the SM with 
enough signihcance, therefore providing strong indirect evidence for physics beyond the 
SM. In either case, the results would greatly improve our understanding of fundamental 
particle physics. In this paper, we perform a study of the constraints on non-universal 
modihcations of the Zbh couplings from prospective precisions at future colliders, 

which, to our best knowledge, is the hrst study of such kind. 


The rest of this paper is organized as follows. In Section we review the current 
constraints on the Zbb couplings and discuss the importance of including the strong cou¬ 
pling constant in the global £t. In Section we compare the precision reaches of the 
three proposed colliders and outline the most important measurements needed for 
improving the Zbh coupling constraints. We then perform a model independent analysis 
to constrain the effective Lagrangian responsible of the modihcations of the Zbh vertex, 
in both the case that the results are SM-like and the one that NP causes a signihcant 


Ref. 


15 


quotes a value of 0.00015 as the theoretical uncertainty for the SM prediction for i?°. Either 


way, the current theoretical uncertainty of is too small to have an impact on the current precision 
data. 
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deviation in the bottom asymmetries. In Section]^ we interpret these constraints on spe- 
cihc NP scenarios including two Higgs doublet models (2HDMs), composite Higgs models 
and the Beautiful Mirror Model. We will compare these constraints to the direct reach 
of the NP particles at the LHC, as well as to the constraints from oblique parameters 
and Higgs precision measurements. In Section we present our conclusions. Finally, in 
Appendix]^ we discuss our treatment of theory uncertainties. 


2 Current constraints on the Zbb couplings 


In this section we present the constraints on the Zbb couplings from the global £t of 
the current precision electroweak data. We follow closely the £t procedure of the Ghtter 
group 14,1^ with the Z-pole data from LEP, SLD 12 along with the measurements 


of the W, top and Higgs masses 24-27 and the hadronic contribution to the running 
hne structure constant 2^ In particular, our procedure uses the input parameters and 
the electroweak and QCD corrections to the electroweak observables of the latest GFitter 
analysis 14 , even though it does not make use of the GFitter public code. 


There are, in fact, a few important differences between our procedure and the one of 
the Ghtter group. First, we include the world average of the strong coupling constant 
as{Mz) (denoted as a 5 (M|)avg.) as a constraint in the htj^ Later in this section, 

we will show that the inclusion of a 5 (M|)avg. has a moderate, but non-negligible, impact 
on the constraints of the Zbb coupling^ Second, the observable sin^ 6*gg(QFB) is not 
included in our £t. sin^0gg(QFB) is a direct measurement of the leptonic effective weak 
mixing angle at LEP, using the charge forward backward asymmetry. Its measurement has 
a strong model dependence and the LEP result explicitly assumes the SM, therefore it is 
difficult to interpret this measurement in the presence of vertex corrections. In practice, 
this observable is not precisely measured and has a small impact in the global £t. To 
summarize, the SM free parameters considered in our £t are Mh , Mz , rutop , and 

as{M‘^). The additional observables included in our £t are F^, Apg, M;, Ac, 


19 


23 


^For a global fit in the context of the Standard Model effective field theory, see e.g. 

®It should be noted that there also exist non-trivial bounds on the Zbb couplings from hadron colliders, 
as pointed out in e.g. Ref. 29 , although the precision can not compete with the one obtained from lepton 
colliders. 


^The value we use is = (2757 ± 10) x 10 

included Q!s(Mf )avg. but did not comment on its impact on the Zbb couplings. 
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Ab, ^fb) ^fb! and Tw- 

The coupling of the Z to left and right handed bottoms, denoted by gib and g^b, are 
given in the following interaction term, 

^ 3 —Z^{gLbbL7^bL + gRbiR'y^bn ), ( 2 . 1 ) 

Cw 

where cw = cos9w, with 9\y the weak mixing angle, and g is the 517(2) gauge coupling. 
Through out this paper we shall use 6gLb and 6gjib to parameterize the modihcation of 
the Zbb couplings, dehned as 


9Lb — + SgLb , gRb — 9^ + Sgub , 


( 2 . 2 ) 


where gf^ and g^ are the SM predictions for g^b and g^b, which at the tree level are 
given by 

=-1/2+ 5^3--0.42, 


SM,0 
9Lb 


9|“'” = s^k/3^ 0.077. 


At the tree level, R^, Ab and Apg at Z-pole can be written as 




9Lb + 9Rb 

Ei9l, + 9%)^ 


(2.3) 


(2.4) 


where denotes a sum over all quarks except the top quark, and 

Q 


Ab = 


9lb 


9Rb 


gib + 9Rb 


A0,b _ 3 - - _ ^ gie 9Re9Lb 9Rb 

^FB A A 2 , 2 2 < 2 

4 4 gfpg + g^ij + gj^f^ 


(2.5) 


These expressions will be modihed once loop corrections are included. 

It should be pointed out that the measurements at Z-pole alone could not determine 
the signs of gib and gRb- The off-peak measurements can resolve the sign ambiguities due 
to interference of the Z diagram with the s-channel photon diagram. However, as pointed 


out in Ref. 32 , the LEP data at scales different from mz are limited in statistics and 
could dehnitely resolve the sign of gib, but not the one of gRb- Nevertheless, NP theories 
able to flip the sign of gRb are typically in tension with other EW precision data such 
as the constraints on S and T parameters 16,32 . Future lepton colliders will collect a 


large amount of data at higher scales, which should completely resolve this ambiguitjj^ 


®As an example, the value of Tpg at around 240 GeV is changed by ~ 0.2 if the sign of gub is 
flipped with respect to the SM prediction. On the other hand, the proposed CEPC run at ^ 240 GeV 
would collect 5 ab~^ of data over ten years with two detectors 1^, resulting in a statistical uncertainty of 
^ 0.0003 for Apg, which is sufficient for resolving the sign of gRt as long as enough events are left after 
the event selection and the systematics are under control. 
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In this paper we do not consider the possibility that gRb has the opposite sign of its SM 
prediction. 

The potential NP that modifies the Zbb couplings can also change other EW observ¬ 
ables. To capture the most relevant corrections without relying too much on the model 
assumptions, we will consider NP scenarios that contribute to the oblique parameters S 
and T along with the modified Zbb couplings. Therefore, our minimal model assumption 
is SM together with S, T, 5gLb and Sg^b treated as free parameters. For later convenience 
we will denote it as {SM+S,T,6gLb,bgRb)- 

With the model assumptions and £t procedure described above, we obtain the con¬ 
straints on dgib and SgRb, shown in Fig. The blue (orange) region corresponds to a 
confidence level (CL) smaller than 68% (95%), while the green dot is the SM prediction 
{SgLb = bgRb = 0). In addition, in the left plot we show the individual constraints from 

(red) and the combination of Ab and Apg (cyan), for which the parameters other than 
SgLb and SgRb are set to the best-fit values. This verifies that and Apg(Mfe) are the 
most relevant measurements for constraining the Zbb couplings. Given that bgit and 6gRb 
are relatively small, the change in Apg and Ab can be expanded in terms of 5gib and 
SgRb- Keeping the first order, while fixing the other parameters to the best-fit values, we 
obtain 


- 0.78 Sgib+ 0.14:6gRb, 

Mp’g ^ - 0.0345gLb-0.185gRb^ 0.11 X 6Ab. (2.6) 

As shown in Table the current experimental (la) uncertainties are 0.00066 for i?® and 
0.0016 (0.020) for Apg {Ab). These numbers together with Eq. (2.6) provide a good 


analytical understanding of the Zbb coupling constraints. First, R^ is more constraining 
than Ap 3 (M;,) and leads to a large positive correlation between Sg^b and SgRb. Second, 
is numerically more sensitive to Sgib while Apg(M6) is more sensitive to SgRb. 

As mentioned previously, in our fit we have included as{Mz)a.vg. (world average) as 
a constraint in the global fit. as{Mz)a.vg. includes several different measurements, but is 


dominated by the lattice calculation 24,30 . More explicitly, to avoid double counting. 


we use the PDG world average excluding electroweak precision test (EWPT) results [24] , 
which is 

as{Mz)s,vg. = 0.1185 ±0.0005 (world average w/o EWPT result) . (2.7) 


The Gfitter group (I4,18 does not include this constraint in the fit, since the global fit 


for the SM and also for the oblique parameters S and T is not sensitive to Q! 5 (M|)avg.. 
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Figure 1: The preferred region in the {Sgiby^gm) plane obtained by the global fit to 
{SM+S,T,SgLb,SgRb) with current data. The blue (orange) region corresponds to the 
68% (95%) CL and the green dot is the SM prediction {5gLh = ^gRb = 0). Left: The 
individual constraints from and from the combination of Ab and are shown in red 
and cyan, respectively. For these curves, the parameters other than 5gib and SgRb are set 
to the best-£t values. Right: The purple contours show the preferred regions for which 
a 5 (M|)avg. is not included in the £t. The solid line corresponds to a 68% CL and the 
dotted line corresponds to a 95% CL. 


One can simply extract the value of asiM^) from the EW global £t, which is in good 
agreement with assuming the SM {+S and T). This can be seen by our hrst 

two results in Table However, the extraction of from EW global £t has some 

model dependence and, with our model assumption, (SM+S', T, Sgib, SgRb), the agreement 
with Q!s(M|)avg. becomes a bit worse (but still < la). This result is shown in the last row 
of Table 1^ This suggests that including a 5 (M|)avg. in the fit can have some impact on 
the Zbb coupling constraints. Indeed, as shown in the right plot of Fig. the inclusion of 
as{Mz)e,vg. has some small but non-negligible effect on bgib and on its correlation with 
SgRb- This is because as{Mz)s.wg. has a stronger effect on than on Hpg and prefers 


smaller 5i?°s, leading to an increase in bgL^ and a smaller decrease in bgRj, (see Eq. (2.6)). 


In the future, while both the precision of the Z-pole data and as{M^)avg. will be improved, 
a 5 (M|)avg. will at least provide an important consistency check and it will be interesting 
to include it in the global £t for the Zbb coupling constraints. 

To summarize this Section, in Table we list both the individual values of and the 
correlations among S, T, Sgib, bgRb obtained from the global £t of the current precision 









model 

as{Mz) from EW global £t 

SM 

0.1185 ±0.0026 

SM+(^,T) 

0.1180 ±0.0027 


0.1153 ±0.0035 


Table 2: Prediction for as{M‘^) from EW global fit with different model assumptions and 
without the world average a 5 (M|)avg. = 0.1185 ± 0.0005 as a constraint. 



S 

T 

^9Lb 

^9Rb 

A 

-0.047 ±0.097 



T 

0.91 

0.015 ±0.077 


^9 Lb 

-0.34 

-0.23 

0.0030 ±0.0015 


Sgm 

-0.40 

-0.30 

0.91 

0.0176 ±0.0063 


Table 3: Best £t values (±1 ct) of and correlations among S', T, dgib, from current pre¬ 
cision electroweak data, with model assumption being (SM-I-S, T, 5gLb, ^gm)- A Gaussian 
distribution is assumed. 

electroweak data, with model assumption being (SM-I-S', T, Sgih, ^gm)- A given NP model 
may contribute to both S', T and Sgib) and to constrain the model one should in 

principle include all four parameters. However, as shown in Table the correlation 
between the two groups, (S', T) and {Sgiby^gm), are not very strong, and we expect a 
similar behavior at future colliders, given that the relative improvements are not extremely 
different for different observables. For simplicity, in the next Section we shall focus on the 
constraint on Sg^b and 6gRb and marginalize over S and T. We refer the reader to other 
literature, e.g. Ref. [^, for prospective constraints on S and T at future colliders. 


3 Zhh coupling constraints from future e+e colliders 


Future e’''e“ colliders will be able to significantly improve the precision of the measure¬ 
ments at the Z-pole thanks to a much larger statistics. The reaches have been estimated 
in the Technical Design Report (TDR) for ILC (^, the TLEP whitepaper for FCC-ee 
and the preliminary Conceptual Design Report (preCDR) for CEPC [^. However, these 
estimations usually either contain only a subset of EW observables, or have combined 
several observables into one {e.g. the effective leptonic mixing angle sin^^g), and are 
therefore not straight forward to apply in our study. In addition, some of the estimations 
are rather preliminary, having strong dependence on the assumptions for systematic un¬ 


certainties and whether or not beam polarization will be implemented. In Section 3.1, we 
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outline the key observables that are needed for improving the Zbb coupling constraints 
and try to estimate their precision reach at the three future colliders. Using these estima¬ 
tions, we proceed to study the constraints on the Zbb couplings by the method of global 


£t, and the results are shown in Section 3.2 and 3.3 


In our study we consider the following benchmark scenarios for the three colliders: 


• CEPC with a relative conservative estimation for the systematic uncertainties and 
with a statistics of only 2 x 10® Z events. While beam polarization could be a 
potential option for the run at the Z-pole, here we assume that it is not implemented. 

• CEPC+, which is CEPC with a more aggressive estimation for the systematic 
uncertainties, and assuming 10^® Z events. 

• ILC, with a lower statistics (10® Z events), but with beam polarization. 

• FCC-ee with 10^^ Z events and beam polarization. 


3.1 Precision of the EWPOs at future e“^e colliders 

The observables directly related to the Zbb couplings are R^, dpg (measured without 
beam polarization) and Ab (measured with beam polarization). However, the three ob¬ 
servables also have explicit dependence on the effective weak mixing angles, sin^ (for 
leptons) and sin^ (for bottom)]^ In particular, Hpg is quite sensitive to sin^ as 
it is proportional to Ae- In fact, at present, the LEP measurement of Hpg provides one 
of the best determination of sin^^g^ assuming SM (the other one being Alr from SLD). 
Therefore, to extract the Zbb couplings from Hpg, it is important to obtain an indepen¬ 
dent determination of sin^ while the most precise such determination is provided by 
the leptonic asymmetry observables. On the other hand, and At are numerically not 
very sensitive to sin^ 

Without beam polarization, the forward-backward leptonic asymmetry Hpg (= \A^) 
can be measured. In addition, a measurement of Ai {Ae and Ar) can be obtained using the 
average hnal-state longitudinal r polarization and its forward-backward asymmetry [I^ , 
which we denote as Ai{V^). With beam polarization, the left-right asymmetry Arr (= Ae) 
can be directly measured, but it is rather irrelevant in terms of the Zbb coupling constraints 


10 


This can be seen from the tree-level expressions in Eq. (2.41 and Eq. (2.5). The effective weak 


mixing angle for leptons (sin^ and bottom (sin^ 9^^) are not exactly the same due to different loop 


contributions. 
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as Ab can be directly measured as well, and we have checked that the impact of the 
improvement of Am on the Zbb coupling is rather negligible. However, can still be 
helpful for constraining the Zbh couplings in a global fit, for example, in the case that 
two (or more) colliders are built and only one of them have beam polarization, similar to 
the situation of LEP and SLC. 


It is worth noting that, apart from a number of additional observables are also 
sensitive to the coupling combination through the dependence on the total 

hadronic decay width, among which R^, which is the ratio of the total hadronic Z decay 
width and the Z decay width to one lepton species, is relatively well measured and provides 
the best sensitivity. We hnd that a signihcant improvement of the precision of can have 
a significant impact on the Zbb coupling constraints. In particular, with the estimated 
precision reach at FCC-ee shown later, the measurement of R^ turns out to be more 
constraining than the one of R^ to the Zbh couplings. However, the constraint from R^ 
depends strongly on the assumption that the coupling of Z to other fermions are SM-like, 
and one should be cautious when applying the Zbh coupling constraints to a model for 


which this assumption is not true. In the end of Section 3^ we will also show the results 
for FCC-ee with a more conservative estimation for the precision of R^. 

To obtain the estimation of the precision reach of the several observables, the following 
procedure is performed. For each observable, we use the estimation in the corresponding 
literature, if it is provided. In particular, if a range of values is provided, we choose the 
more conservative one. If the estimation is not provided in the literature, we estimate 
the precision with the following strategies: we assume that the systematic uncertainties 
at CEPC is a factor of 1/3 the ones at LEP and the systematic uncertainties for the 
scenario “CEPC-l-” is reduced by an additional factor of 1/2 from CEPCp^ For CEPC 
and CEPC-I-, we assume the statistical uncertainty simply scales with 1/\/N, where N 
is the total number of Z events expected to be collected. Additionally, for ILC, Ref. 
does not provide an estimation for the uncertainty of i?/ (Ai?/), for which we adopt the 
estimation in Ref. 18 by Gfitter, AR/ = 0.004. Finally, for the FCC-ee, Ref. does 
not provide an estimation for the uncertainty of Ab (A^;,). We naively scale it from the 


estimation for the ILC, assuming 
at FCC-ee. 


I ILC 


AAr.p 


IFCC 


_ee, which gives A^^ 


0.00021 


exception of this is the CEPC systematic uncertainty for R10(0.006), which we deduced from 
the total uncertainty (0.007) in the preCDR and then scale to obtain the systematic uncertainty at 


CEPC-L (0.003). 
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Precision 

Observable 

Current 

CEPC 

CEPC+ 

ILC 

ECC-ee 

n'i 

[0.216291 

0.00066 

(0.00050) 

0.00017 

(0.00016) 

0.00008 

(0.00008) 

0.00014 

0.00006 

(0.00006) 

Ifl 

[20.767] 

0.025 

(0.007) 

0.007 

(0.006) 

0.003 

(0.003) 

0.004 

0.001 

(0.001) 

Ab 

^FB 

[0.0992] 

0.0016 

(0.0007) 

0.00015 

(0.00014) 

0.00007 

(0.00007) 



^FB 

[0.0171] 

0.0010 

(0.0003) 

0.00014 

(0.00010) 

0.00007 

(0.00005) 



MK) 

[0.1465] 

0.0033 

(0.0015) 

0.0006 

(0.0005) 

0.0003 

(0.0003) 



Ab 

[0.923] 

0.020 



0.001 

0.00021 

(0.00015) 

Alr 

[0.1514] 

0.0022 

(0.0011) 



0.0001 

(0.0001) 

0.000021 

(0.000015) 

# of Zs 

~ 2 X 10 '’ 

~ 2 X lO’^ 

-IQi^ 

~ 10^ 

~ 10^^ 


Table 4: The estimated precision reach for the observables (which current experimental 
measurement is shown in the hrst column) most relevant to constrain the Zbb coupling at 
future colliders. The second column shows the uncertainty of the present measurements 
from LEP and SLC, while the other columns show the estimations of the precision reach 
for different future colliders and scenarios. The numbers highlighted in blue are our own 
estimation. In each entry, the number at the top shows the total uncertainty while the 
number at the bottom (in parenthesis) shows the corresponding systematic uncertainty. 
A blank entry denotes an observable that is either not measured or not important for our 
global £t. The last row shows the expected number of Z events that will be collected. 


The estimations for the observables mentioned above are summarized in Table HI A 
similar method is used to estimate the precision reach for the additional EW observables 
not listed in Table even if we have checked that they have a much smaller impact on the 
Zbb coupling constraints. In Table the numbers highlighted in blue are our own new 
estimations for those observables not reported in the literature. In each entry, the number 
at the top shows the total uncertainty while the number at the bottom (in parenthesis) 
shows the corresponding estimated systematic uncertaintj-For comparison, in the 


second column of the table, we show the current uncertainties, taken from Ref. 12 . A few 


comments on the table are in order: the systematic uncertainties mostly dominate except 


for Ab and Am at FCC-ee, and for Apg at CEPC. Without beam polarization, CEPC 


^^For some of the observables, only the total experimental uncertainties are quoted in the literature 
and the entries for the systematic uncertainties are left blank, accordingly. 
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can not measnre Ab and Alr. For this reason, the corresponding entries are left blank. 
ILC and FCC-ee can measnre Ape, v4pg and Ai{V^), bnt the corresponding observable 
with beam polarization, Ab and Arr, can be measnred signihcantly more precisely, so, 
for simplicity, we do not inclnde the former observables. 

A potential issne for the interpretation of fntnre measnrements is the effects of the 
theoretical nncertainties, which could become important if they are much larger than 
the experimental uncertainties. In our study, we assume that the electroweak three-loop 
corrections will be computed in the future. In that case, the effects from the theoretical 
uncertainties on the Zhh coupling constraints are very small and can be safely neglected 
even for FCC-ee. This is either because the theoretical uncertainty is numerically small 
(such as (fth-Rb which is estimated to be a few times 10“^ [^) or the observable itself has 
little effect on the Zhh coupling constraints, such as the top quark mass. The theoretical 
uncertainty of sin^ 9^^ also has little impact, since Ab is not very sensitive to it. More 
details on the treatment of the theoretical uncertainties can be found in Appendix 

3.2 SM-like measurements and constraints on NP 

In this Section, we assume the future experimental results agree perfectly with the SM 
predictions and the estimated precision of future measurements as described in the pre¬ 
vious Section. The preferred regions in the {Sg^b, ^Qm) plane obtained by our global £t 
are shown in Fig. The plot in the right panel is a zoomed-in version of the one in the 
left panel. 

From Fig. it is clear that the constraints on the Z66 coupling are signihcantly im¬ 
proved at the future colliders, even for the relatively conservative CEPC estimation 
(cyan contours), compared to the results of the current precisions (purple contours). With 
beam polarization, ILC (red contours in the hgure) and FCC-ee (black contours in the 
hgure) have better measurements of Ab, which gives a better constraint on 5gRb and also 
reduce the correlation between hgib and SgRb- 

We report the la uncertainties of Sgib and 6gRb as well as their correlation (p) in 
Table 1^ Due to the strong correlation between dgib and 5gRb (in particular at CEPC), 
one need to be careful when using these results to constrain NP models, since in some 
models only one of Sg^b and SgRb can receive a sizable contribution while the other one 
is close to zero. Therefore, in Table we also show the la uncertainties for hgibiSgRb), 
while SgRb{6gLb) is hxed to zero. 
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Figure 2: Preferred regions in the {Sgib, ^gm) plane, assuming SM central values for all 
measurements. The model assumption is {SM+S,T,6gLi,,SgRb), with S,T,6gLb,SgRb all 
treated as free input parameters. The solid and dotted lines are 68% and 95% CLs, 
respectively. The purple contours assume current precision for all measurements. The 
cyan, blue, red and black contours correspond to the estimated precisions for CEPC, 
CEPC+, ILC and FCC-ee, respectively. The plot in the right panel is a zoomed-in 
version of the one in the left panel. 


For CEPC, the estimation for in the preCDR seems to be very conservative, 
suggesting little improvement of its systematic uncertainty from LEP to CEPC (see Table 
1^. A scaling of a factor of 1/3 on the systematic uncertainty would give a value of 
Ai?/ = 0.003 for the total uncertainty. The results for this scenario are shown in Table 
which exhibits a slight improvement. In Tablewe also show the results for FCC-ee with 
a more conservative estimation of Ai?/, also using Ai?/ = 0.003 (instead of Ai?/ = 0.001). 


3.3 Discovering NP through Apg (A*) 

A more interesting possibility is that the long standing Apg discrepancy does come from 
NP, in which case the precision reach at any of the three future e’''e“ colliders should be 
able to rule out the SM with very high signihcance and therefore provide strong indirect 
evidence for physics beyond SM. To illustrate this point, we consider the following two 
scenarios. Scenario I: we assume that the true values for Sgib and 6gRb (denoted by 
Sg^fj and are given by the best £t values of the current data, Sg^f^ = 0.0030 and 
6g^b = 0.0176 (see Table [^. Scenario 11: we assume that the the true values of 6g^f^ and 
Sg^fj are closer to zero, while still being consistent with the current measurements within 
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SgLb 

Sgm 

P 

Sgib {Sgm — 0) 

Sgm {SgLb — 0) 

current 

0.0015 

0.0079 

0.91 

0.00061 

0.0032 

CEPC 

0.00031 

0.0016 

0.87 

0.00015 

0.00079 

CEPC+ 

0.00015 

0.00078 

0.88 

0.000072 

0.00037 

ILC 

0.00016 

0.00059 

0.61 

0.00012 

0.00047 

FCC-ee 

0.000044 

0.00012 

0.42 

0.000040 

0.00011 


Table 5: A comparison of precision reach at different future colliders. The 2nd(3rd) col¬ 
umn shows the la uncertainties of SgLb{SgRb) while marginalizing over 6gRb{SgLb)- The 
4th column shows the correlation (p) between SgLb and Sgjib- The 5th(6th) column shows 
the la uncertainty of Sgibi^gnb) with SgRbi^gib) fixed at zero. We assume future mea¬ 
surements to be in perfect agreement with the SM predictions. A Gaussian distribution 
is assumed. 


AR'^ = 0.003 

SgLb 

Sgm 

P 

Sgib {Sgm — 0) 

Sgm {SgLb — 0) 

CEPC 

0.00031 

0.0015 

0.93 

0.00011 

0.00057 

FCC-ee 

0.000051 

0.00012 

0.32 

0.000048 

0.00012 


Table 6: Same as Table but for CEPC and FCC-ee both with = 0.003, which 
serves as a reasonably optimistic estimation for CEPC and a conservative one for FCC-ee. 

68%CL. As a benchmark point, we choose 6g^fj = 0.0009 and Sg^f^ = 0.0075. In principle, 
one would expect the NP to have non-zero contributions to the S and T parameters, as 
well. We find that changing the central values of S and T of the hypothetical measurement 
within the current constraints has very small impact on the Zbb couplings. For simplicity 
we assume that the hypothetical data agrees with SM other than the modification to gib 
and gnb. 

The preferred regions in the {SgLb, Sgnb) plane are shown in Fig. The two plots 
correspond to the two scenarios described above, and each shows the 99.9999% CL (cor¬ 
responding to 5a for a one-dimensional Gaussian distribution) constraints from different 
colliders. From the figure, it is clear that the SM prediction at zero (denoted by a green 
dot) can be ruled out at 99.9999% CL by all the e+e“ colliders we discuss, even if we 
assume that the future measurements will point towards smaller values of Sgib and Sg^b 
within 68% CL of the current measurements. 


4 Implication on NP models 

In this Section, we analyze the implications of the future measurements of the Zbb cou¬ 
plings on specific NP models. We start with a brief discussion of the constraints on 
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^dRb ^dRb 

Figure 3: The preferred regions in the {Sgib, ^Qm) plane, given by the global £t of the 
future measurements at CEPC (in cyan), CEPC+ (in blue), ILC (in red) and FCC-ee 
(in black). The solid and dotted purple contours correspond to the 68% and 95% CL 
constraints from the current measurements. The two panels correspond to Scenario I 
and Scenario //presented in the text, and each plot shows the 99.9999% CL constraints 
from different colliders with dashed contours. The green dot is the SM prediction {dgib = 
Sgm = 0). 


effective Lagrangians. At dimension 6, the only operators that modihes directly the Zbb 


couplings are (see e.g. 17,21 


Om = i(H^D„H)(bRYbR), 

,j‘D^H)(QYrQ). 




0‘„Q = i{H' 


(4.1) 

(4.2) 

(4.3) 


After electroweak symmetry breaking, these operators lead to a shift in the Zbb couplings: 


^9 Lb — 




^9Rb — 


amv 


(4.4) 


where aubi o,hqj coefficients of the Om, ^hq operators, respectively 

and V is the vacuum expectation value of the Higgs {v = 246 GeV). In Table we present 
the constraints on these operators at the several future machines, assuming that 

ajjQ = a\jQ = am = 1/A^. Scales as large as (20 — 30) TeV can be probed by the future 
measurement of the Zbb couplings. 


Next, we pass to the analysis of specific NP frameworks that can generate some of the 
operators forementioned, including two Higgs doublet models, composite Higgs models 
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current 

CEPC 

CEPC+ 

ILC 

FCC-ee 

A(TeV) 

6.8 

13 

20 

15 

27 


Table 7: 95%CL bounds on the cutoff scale A from different future colliders, using the 
results from Section 


3.2 


and assuming that a^^jQ = = aub = 1/A^. 


and the Beautiful Mirror Model. It should be noted that for natural SUSY with mini¬ 
mal ingredients, the Zbibi coupling receives corrections from loops involving stops and 
Higgsinos, while they are rather small for the future measurements of the Zbb to have an 
impact, as shown e.g. in Ref. [^. 


4.1 Two Higgs doublet models 

As shown in e.g. [^, models with an extended Higgs sector can predict sizable NP 
contributions to the Zbb vertex. In particular, focusing on two Higgs doublet models 
(2HDMs), based on discrete symmetries to avoid flavor changing neutral currents (FCNCs) 
at the tree level, the most important contribution generically comes at the one loop level, 
from the charged Higgs exchange. The sign of the charged Higgs NP contribution to Sgib 
{^9Rb) is fixed and is always positive (negative). In a Type H 2HDM, the contribution to 
SgLb {bgRb) increases at small (large) values of tan/3, since the coupling H^bitR {H^tibR) 
leading to a non-zero 5gLb {bgnb) is proportional to mt/ta.n/3 (mfetan/3). In a Type I 
2HDM, instead, both Sgib and 6gRb increase at large values of tan/Q leading always to 
a NP contribution bg^b ^ In Fig. we show the constraints on the mH± — tan/3 

plane, using the present measurement of the Zbb coupling (in purple) as well as the 
expected more precise measurement at CEPC (in cyan), CEPC-I- (in blue), ILC (in red) 
and FCC-ee (in black). For the figure, we have assumed that the future measurements 
perfectly agree with the SM predictions and we have marginalized over the values of the 
S and T parameters. 

In Type H models, interesting constraints arise at low values of tan (3 for which bgLb 3> 
Type I models, instead, are only allowed in the region with small tan /3 unless 
mH± is very large. If we specify the full spectrum of a 2HDM, including the masses of the 

^^Here we use the tan/3 convention such that the two charged Higgs couplings are proportional to 
mt tan /3 and mb tan /3. 

^"‘in a Type II 2HDM, only very large values of tan/3 can be excluded by the measurement of Sgnb'- 
tan/3 ~ 0(50) (30) for mH± ^ 200 GeV at CEPC-I- (FCC-ee). At large tan/3, the exact exclusion 
depends also on the neutral Higgs spectrum, since in this case the neutral and charged Higgses give a 
parametrically similar contribution to SgLb and Sgub- 
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Type II 



Type I 



Figure 4: Constraints on the charged Higgs parameter space in a Type II 2HDM (left 
panel) and in a Type I 2HDM (right panel). In purple is the constraint we obtain using 
the present uncertainties on the EWPOs; in cyan, blue, red and black are the constraints 
expected with the future measurements at CEPC, CEPC+, ILC and FCC-ee, respectively. 
In the Type II 2HDM, the region below the curves is excluded. In the Type I 2HDM, the 
exclusion is above the curves. We assume that the future measurements perfectly agree 
with the SM predictions and marginalized over the values of the S and T parameters. 


neutral scalar H and pseudoscalar A, as well as the mixing angle a — (3 between the two 
doublets, S and T are not free parameters. In general, the constraints will be stronger. 
As an example, in a type II 2HDM, if we £x tua = 'mu = 200 GeV, mH± = 150 GeV and 
cos(a — (3) = 0.1, then the model gives S ~ 0.02 and T ~ 0.04. For this specihc set of 
parameters, the GEPG bound on tan/3 is ~ 2, instead of ~ 1.8, as presented in the left 
panel of Fig. (see the cyan dashed line). 


Presently, LHG charged Higgs searches almost totally exclude charged Higgs bosons 

There are no LHG searches 


with a mass below the top mass in Type H 2HDMs 33 


at around the top mass, for 160 GeV < m[j± < 180 GeV up to date. Above 180 GeV, 
constraints are rather weak and cover only models with large values of tan (3 (tan (3 > 40), 
for which the production cross section of the charged Higgs in association with a top 
quark is in the 0{1) pb range. In this regime, the two most important bounds come from 
the searches for —)■ tu 33 and for th 34 . At the 14 TeV LHG, also charged 


Higgs boson with mass above the top mass will be relatively well probed. In particular, 
searches for th will have the potential to probe tan/3 < 3 and tan/3 > 15 for 

mH± ~ 500 GeV at the High Luminosity (HL)-LHG 35 . Gomparing to our results of 
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Fig. I (left panel), we see that constraints from fntnre measnrements of EWPOs can be 
complementary to direct searches for Type II 2HDMs, being able to probe low valnes of 
tan/3 even for mH± 3> m*, as well as the challenging region 160 GeV < mH± < 180 GeV, 
presently not covered by direct searches. 

Finally, one can interpret the searches for charged Higgs bosons in terms of Type I 
2HDMs. Below the top mass, only a small region with tan /3 < 1 has not been yet probed 
by the —>■ tv search. Above the top mass, the exclnsion is very week and is not 

covering any part of the plane shown in Fig. (right panel). At the HL-LHG, this region 
will be very well probed by a tb search, with potential exclusions for the entire 

range of tan (5 presented in the hgure, up to mH± ~ 500 GeV. 


4.2 Composite Higgs models 


Gomposite Higgs models usually predicts a large correction to the Zbibi coupling, since 
a sizable mixing between the third generation quarks and the strong dynamics is needed 
to generate the large top mass. The correction to the Zbpjbn coupling is usually much 
smaller, unless one specifically extend the fermion sector to generate a large correction 
{e.g., as in Ref. 1^). It was pointed out in Ref. 10 that an 0(4) symmetry, which 


is the SU{2)l ® SU{2)ji symmetry, analogous to the custodial symmetry protecting the 
weak isospin, with the addition of a left-right parity Plr, could be used to protect the 
Zbibi coupling, such that a natural composite Higgs model can be consistent with EW 
precision constraints. Nevertheless, in realistic models there still exist corrections to the 
Zbibi coupling because 1) the Plr symmetry can only protect the Zbibi coupling at zero 
momentum and 2) there are also several contributions that explicitly break Prr- These 
corrections could become relevant if the constraint on the Zbibi coupling is significantly 


improved at future e^e colliders. Ref. 37 estimates the size of different contributions 


to the Zbibi coupling in minimal composite Higgs models with custodial protection. 


(Also see Ref. 38 for a recent review.) While the magnitudes and signs of different 
contributions are rather model dependent, the leading correction usually comes from Plr 
breaking effects of fermion loops and is 


^gib 


Vt 




( 4 . 5 ) 


IGvr^ 

where yt is the top Yukawa coupling, rup is the mass of the p meson which cuts off the 
loop correction, and rri/^ is the mass of the 4-plet composite quarks which is essentially the 


mass of the (lightest) top partner up to some corrections from mixing. Taking Eq. (4.5) 


19 









Figure 5: Constraints in the {gp/giij,f) plane assuming the only correction to the Zbb 
couplings are given by Eq. (4.5). Each contour represents the 95% CL constraint and the 
region in the top-left side of the contour is excluded. The grey horizontal line corresponds 
to gplgij, = 3/0.7 . In purple is the constraint we obtain using the present uncertainties 
on the EWPOs; in cyan, blue, red and black are the constraints expected with the future 
measurements at CEPC, CEPC-I-, ILC and FCC-ee, respectively. 


as an equality, the results in Table (assuming 5gRh = 0) can be interpreted in terms 
of constraints in the {gp/gip^f) plane where gp = rrip/f and g^ = m^/f. This is shown 
in Fig. where each contour represents the 95% CL constraint and the region in the 
top-left side of the contour is excluded. The grey horizontal line has gpjg^ = 3/0.7 which 
corresponds to the benchmark point nip = 3 TeV and = 700 GeV of Ref. . 

Since gp/g^ is typically bounded to be a few times one, future e+e“ colliders can 
constrain / to be at least a few TeVs thanks to the measurement of the Zbb couplings. 
This is comparable to the constraints from the direct searches of top partners at the 
next run of the LHC, given that the mass of the top partner can not be much larger 
than / in order to obtain the correct Higgs mass 39 . The constraints from the Zbb 


couplings is signihcantly stronger than the ones from oblique parameters but weaker than 
the ones from Higgs precision measurement, and in particular from the HZZ vertex, 
quoted in Ref. [^. The latter can, in fact, constrain / at the level of ~ 2.8 TeV (CEPC) 
and ~ 3.9 TeV (FCC-ee) at 95%CL. Other studies of future constraints on composite 
Higgs models can be found in Ref. 40,41 . To conclude, while the constraints from Zbb 
couplings has a stronger model dependence, it is complementary to the constraints from 
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oblique corrections, Higgs precision measurements and direct searches and can be very 
helpful for the discrimination of different composite Higgs models. 


4.3 The Beautiful Mirror Model 


In the Beautiful Mirror Model proposed in Ref. 32 , the modihcations to the Zbb couplings 


are caused by the mixing of the bottom quark and new heavy vector-like quark (s)p^ In 
the simplest case with only one bottom-partner B, the shifts in the couplings are given 
by ^9Lb = (ts + l/2)s| and Sg^b = where is the charge of the new bottom quark, 
B, under SU{2)l and sl(r) is the sine of the left(right)-handed mixing angle between 
the SM b quark and B. To obtain a shift in g^b, the new quark B can not be a SU{2)l 
singlet. If B has 6gRb would be negative and one would need a large shift in 

the coupling to flip its sign, g^b ~ order to resolve the Hpg discrepancy. Such a 

large shift requires a very light B quark and a large custodial symmetry breaking and has 


been almost completely probed by LHC direct searches for vector-like quarks 16,43,44 


A more appealing choice is fs = that can lead to a good £t of the EWPOs, without 
a too light B quark, thanks to a positive contribution to g^b- This scenario was denoted 


as the “Top-less Mirror” in the original paper 32 , since there is no top-partner in the 


new doublet quark. The new doublet quark alone can not simultaneously generate a 
sizable enough SgLb, but this can be easily achieved by introducing an additional singlet. 
Therefore, we discuss an extension of the SM with the following vector-like quarks. 


'^L,R — 


B 


L,R 


^(3, 2,-5/6), 
(3,1,-1/3), 


(4.6) 

( 4 . 7 ) 


where the numbers in the bracket denotes representations under SU{3)c, SU{2 )l, and the 
U{1)y hypercharge. The relevant terms in the Lagrangian are given by 

— C D + M 2 BlBr + UiQLHbR + hlQrHBr -|- yR^iBbR + h.c., (4.8) 

which leads to the following 3x3 mass matrix AAr for the bottom-like quarks while the 


Mr — 


X is 

simply given 

by Ml, 


0 



[yr 

Ml 

0 , 

Mx 

lo 

0 

M2 



(4,9) 


^®See also 
quarks. 
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for a more recent analysis of the fit of the Zhb couplings in models with vector-like 
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where = ^- The shifts in the Zbb couplings are thus given by 


SgLb = 


XL 

‘IMl ’ 


— 


XL 

2Ml 


(4.10) 


The new quarks contributes also to the T parameter through fermion loops. Ignoring the 
small bottom mass and the higher order terms, the T parameter is given by 

MimX 


T 




^ ^ghb^i ^^gLb-j^2 


m 


log 


top 




m 


top 


(4,11) 


If SgLb and 6gRb are fixed to non-zero values, as preferred by the current EW measure¬ 
ments, the T parameter will increase if one increases the mass scales of the new quarks, 
as shown in the left plot of Fig. This is because larger Yukawa couplings are needed if 
one wants to fix bgib and Sgub while increasing Mi and M 2 , as shown in Eq. (4.10). Since 
the desired value of Sgjib is much larger than the one of bgib^ fhe T parameter is more 
sensitive to Mi. The contribution to the S parameter is very small in this model, and we 
assume it to be zero for simplicity. The right plot of Fig. (solid blue line) shows the value 
of (measured from the minimum that can be reached in this model) as a function 
of Ml, obtained from the global fit of the current EW precision data, while other model 
parameters are chosen to minimize Ay^ for a given value of Mi. At small Mi (< 1.4 TeV), 
one could tune the value of M 2 to obtain the best-fit value of the T parameter 0.05 
if fixing S' = 0). For larger values of Mi, the agreement with data is significantly worse 
due to the tension between the T parameter and Sg^b- For comparison, we also show 
the Ay^ curve for CEPC (dashed red curve), assuming that the central values of S, T, 
SgLb, Sgnb are the same as the current ones: CEPC will tightly constrains Mi to be below 
~ 1.6 TeV. 

An additional upper bound on Mi can be obtained by requiring the theory to be 
perturbative. In particular, Eq. (4.10) implies yn ~ 2y/SgRbX- Assuming the true values 
oi SgLb and SgRb are within 
and hence 


CL of the current measurements, we have SgRb > 0.0075 

Ml 


Ur 


> . (4.12) 

~ 1.4 TeV ^ ^ 

Therefore, Mi can not be much larger than a few TeV without the Yukawa coupling, 

violating perturbativity bounds. 


The LHC is directly searching for the mirror quarks. As pointed out in Ref. 44 , the 
charge —4/3 exotic quark X decays to b and W with the same sign of electric charges, 
which is extraordinary in theory but hard to capture in experiments, since it is very hard 
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Figure 6: Left: The value of T parameter in the (Mi, M 2 ) plane, while 5gLb and 5gRb are 
hxed to the current best £t values, Sgib = 0.0030 and dgm, = 0.0176. Right: as a 

function of Mi, obtained from the global fit of the EW precision data to the model while 
marginalizing over other model parameters. The Ax^ is measured from the minimum 
X^ that can be reached in this model. The solid blue line is obtained from the current 
data; the dashed red line is from hypothetical CEPC data, while we naively assume that 
the central values of S', T, 6gn,i are the same as the current ones. The horizontal 
axis starts at 912 GeV, the current 95% CL bound from the direct searches of vector-like 
quarks at the 8 TeV LHC. 


(if not impossible) to measure the charge of 6-jets. As such, the strongest bounds on X 
come from searches of top partners decaying to b W. The recent CMS analysis 45 sets 
a lower limit of 912 GeV at the 95% CL for a pair produced top vector-like quark with 
100%BR to bW, using 8 TeV data. Currently, this is the most stringent constraint on 
Ml. There also exist bounds from bottom partner searches {e.g. Ref. |^) which are 
slightly weaker. 

The current constraint from the LHC is consistent with the one from EW precision 
data shown in the right plot of Fig. The situation may get much more interesting in the 
future: the bounds on the mirror quark masses are expected to reach (2-2.5) TeV at the 

using the single production channel. The HL-LHC with 3000 fb~^ 


14 TeV LHC 


44 


47 


data could probably push the bound further to above 3 TeV, if the data is consistent 
Such a bound would be in tension with the current EW precision data. 


with SM 48 


If no mirror quark is found during the LHC run, it could be an indication that either 1) 
the LEP Apg discrepancy is not due to NP or it comes from some NP other than the 
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Beautiful Mirror Mode or 2) the underlying NP is some extension or modification of 


the “minimal” Beautiful Mirror Model which evades the constraints from direct searches, 
still producing a good £t to EW precision data. 

Finally, we comment on the fact that the Beautiful Mirror Model predicts a signihcant 


modihcation to the hbb and hgg couplings 16 
Vhbb 


iVhbb) 


SM 


- 1 - ^Sgib + Sgm) , 


9hgg 


idhgg) 


SM 


~ 1 + 2{5gLb + SgRb) ■ 


(4.13) 


The current best-£t values of 6gLb and 6gRb predict a ~ 4% deviation on the hbb and hgg 
couplings which will be detectable at a future Higgs Factory [T]-[^. 


5 Conclusion 


Precision measurements of SM couplings are the central goal of future lepton colliders. 
Such measurements are complementary to direct searches at high energy proton colliders. 
In this paper, we have extracted the constraints on possible modihcations of the Zbb 
couplings from the SM predictions by performing global hts of both the current precision 
data and the prospective data from future e+e“ colliders. We pointed out that the world 
average value of the strong coupling constant contains non-trivial information and should 
be included in the global £t of models with non-zero NP contributions to the Zbb vertex, 
which has not been pointed out elsewhere. For the future colliders, we summarized the set 
of observables most important for improving the Zbb coupling constraints and compared 
the precision reaches at CEPC, ILC and FCC-ee. We studied both the case that the 
results are SM-like and the one that the Zbb couplings deviate significantly from the SM 
prediction as suggested by the LEP Hpg discrepancy. For the latter case, we showed that 
even if we assume that the future measurements will point towards smaller values of 6gib 
and 6gRb within 68% CL of the current measurements, any one of the proposed e+e“ 
colliders will be able to rule out the SM with more than 99.9999% CL, equivalent to 5 
standard deviations for a one-dimensional Gaussian distribution. 


Finally, we studied the implication on NP models from the improvements of the Zbb 
coupling constraints. We hrst considered generic 2HDMs, in which the limits from pre¬ 
cision Zbh measurements are complementary to those from the LHC searches of charged 


Another possible solution would be a Z' near Z-pole, as proposed in Ref. 


49 


50 


col¬ 
liders could perform a much better scan around the Z-pole compared with LEP and provide significantly 
better discriminating power between the SM and this scenario. 
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Higgs bosons. In particnlar, fntnre measnrements of the Zbb vertex will be able to 
test Type II 2HDMs at low valnes of tan/? (tan/? < 1 — 4), even in the mass range 
160 GeV < mH± < 180 GeV, presently not covered by direct LHG searches. We then 
considered composite Higgs models, where deviations in the Zbh conplings are generically 
expected. Measnrement at fntnre lepton colliders can probe composite resonances with 
masses of mnltiple TeVs, possibly beyond the reach of the LHG. Finally, in the literatnre 
there have been new physics models motivated by the long standing LEP anomaly in 
Hpg. As an example, we considered the so called Beantifnl Mirror Model in which new 
fermions mix with the SM bottom qnark. We hnd that, for the minimal scenario, fntnre 
precision measnrements of Zbb conld pnt a strong npper bonnd on the new fermion masses 
at aronnd 1.6 TeV, if the measnred central valnes of S, T, bgib and 6gRb are the same as 
the cnrrent ones. Fntnre LHG searches of vector-like qnarks will be able to fnlly probe 
this mass range. 

We have also shown that the particnlar models considered in this paper generically 
predict deviations in the Higgs conplings, which can also be measnred very precisely 
at the Higgs factory stage of the lepton colliders. The interplay between Higgs and 
Z precision measnrements will be very valnable at fntnre machines, in extracting 
maximal information abont new physics. 
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A Theoretical uncertainties 

We follow Ref. for the treatment of theoretical nncertainties, which we discuss in this 
Appendix. Given a set of model parameters 6 , the true value of a particular observable 
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X is predicted to be some certain value, Xth, up to some uncertainty, 6. This uncertainty 
could come from the omission of higher order terms in a hxed order calculation, or the 
subtlety in the dehnition of certain quantities {e.g. the top mass). As such, it is strictly 
speaking not a statistical quantity and there is no good reason to assume it follows a 
Gaussian distribution. We assume x takes a probability density function h{x; 6) that is 
flat within Xth ± 5 and zero elsewhere. 



if 

|x - XthI 

[o 

if 

|x - XthI 


J • (A.1) 

On the other hand, we assume the measured value, x^ea, takes a Gaussian distribution 
5'(3^mea|a;) Centered around the true value x with standard deviation a. 

The distribution of the measure value Xmea given a set of model parameters 6 is thus 
obtained by convolution: 


f{Xme 3 .-, 0 ) = jdx g{x^eJ\x)h{x-,6) 

nea — Xth + d 




^2, 


a 


erf 






^2, 
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(A.2) 


which reduces to the Gaussian distribution in the limit 5 —?• 0. Eq. (A.2) can be im¬ 


plemented in a global £t with N observables with a modified function (assuming no 
correlation) 


N 

Xmod 


j=l L 


2 log ( — ( erf 


Mi - Oi + 6, 


a/2ct,- 


Mi - Oi 


erf ( -2 j j j - 2 log 

(A.3) 

where for each observable j, Mj is the measured value, Oj is the predicted value, aj 
is the experimental uncertainty and 6j is the theoretical uncertainty. and Apg(^b) 
are directly related to the Zbb couplings, and their theoretical uncertainties are most 
important to our study. The theoretical uncertainty of {6thRb) is estimated to be 1.5 x 
10““^ from two-loop diagrams without closed fermion loops and higher-order contributions 
15 and could be reduced to a few times 10“® assuming the 2-loop and 3-loop computations 


will be completed in the future [^. Naively, one would expect it to have an impact, 
especially for FGG-ee which will be able to measure R^ to a precision of about 6 x 10“^. 
However, even with a conservative estimation, dthRb = 5 X 10 and for FGG-ee, we find 
the change of the total uncertainty from the inclusion of the theoretical uncertainty rather 
small, as the 68% GL bound changes from 6 x 10“® to 6.65 x 10“®. The corresponding 
probability density functions are shown in Fig. for which we have set the central value 
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to zero and scaled up the uncertainties by 10® for convenience. Given that the estimations 
of the future experimental uncertainties are still very preliminary, we ignore this small 
effect due to the theoretical uncertainty of in our study. 


Probability Density Function 



Figure 7: Probability Density Functions (p.d.f.) for theoretical uncertainty (5th = 5, ex¬ 
perimental uncertainty (Jexp = 6 and the combined total uncertainty, assuming a zero 
central value. The blue solid line shows the p.d.f. with the experimental uncertainty only, 
which follows a Gaussian distribution. The red dashed line shows the p.d.f. with exper¬ 
imental and theoretical uncertainties combined with Eq. (A.2). The light vertical lines 
shows the corresponding 


and 95% GL bounds. Without the theoretical uncertainty, 
the 68% and 95% GL bounds are ±6 and ±12; with the theoretical uncertainty, the 68% 
and 95% GL bounds are ±6.65 and ±13.0. 


With (longitudinal) beam polarization. At can be directly measured. Without beam 
polarization, Apg can be measured, which is related to Ab by Apg = ^AeAb- The value of 
Ae can be extracted experimentally by the measurement of Apg. Therefore, the theoretical 
uncertainty of App also only comes from At, which can be parameterized by the theoretical 
uncertainty of sin^ 6*®^ (The overall form factors of gLb and g^b cancel in the ratio.) Ab is 
numerically not very sensitive to sin^ 0®^. At the leading order in the SM, one has 

SthAb ^ -0.64 (5th sin^ 9^^ , (5thApB -0.070 (5th sin^ 9^^ , (A.4) 


where <5th denotes the theoretical uncertainty of the corresponding quantity. Even with the 
current theoretical uncertainty of sin^ 6*®^, which is about 5 x 10“® 51 , 6thAb and (5thApg 


are much smaller than the future experimental precisions and can be safely ignored. 


The effects of the theoretical uncertainties of other quantities, such as the top mass 
and W mass, are important in general {e.g. for the S and T parameters, as pointed out in 
Ref. [^) but do not directly affect the Zbb couplings. We implemented these theoretical 
uncertainties and found that the changes of the Zbb coupling constraints are negligible. 
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